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§ . Abstract 

A computer algebra algoritm for solving the quantum Yang-Baxter equation is pre- 
£Nj ■ sented. It is based on the Taylor expansion of i?-matrix which is developed up to the 

order A 6 . As an example the classification of 4 x 4 i?-matrices is given. 
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By the Yang-Baxter equation (in the braid group form) is called the following equation []IJ: 
Rn(X - fi)R 2 sWRM = iW)i?i 2 (A)i? 23 (A - n), (1) 

on the N 2 x A^ 2 matrix -R(A) called i?-matrix. Here for an arbitrary N 2 x A^ 2 matrix X the 
iV 3 x A^ 3 matrices and X23 are denned as X12 = X ® In, X23 = In <£> X, where In is a 
unit N x N matrix. 

In this paper we study only the regular solutions of the Eq. (1). These solutions play 
an important role in the Quantum Inverse Scattering Method and additionally satisfy the 
following condition: 

X\ R(0) = C-I N2 , (2) 

5-j 

where C is a numerical constant. 

The system (1),(2) is invariant under multiplication of R{\) on the arbitrary regular func- 
tion /(A) satisying the initial condition /(0) = 1: 

J2(A)-/(A)J2(A). (3) 

More complicated are the SL(N,C) symmetry ||T|: 

R(X) -> M- 1 ® M~ 1 R(X)M <g> M, M e SL(N,C), (4) 
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and two discrete symmetries: 

R(X)^R t (X), (5) 
R(X) -> PR(X)P. (6) 

Here P is the iV 2 x iV 2 exchange matrix defined as operator in <g> C N by the formula: 

■fy.fcj — 8il$jk- (7) 

Our method of solving the Eq. (1) is based on the Taylor expansion of i?-matrix: 

oo 1 

R(X) = ^- J R (n) A n . (8) 

n=0 ' 

If we take for simplicity in the Eq. (2) C = 1 then it follows that R^ = I. (From now we 
shall write I instead of In 2 )- The matrix R^ plays a key role in our approach. As the local 
Hamiltonian density of the integrable model corresponding to R it is usually denoted by H. 
Taking in mind applications to the Inverse Scattering Method we shall study only nontrivial 
Hamiltonian densities. These ones can not be represented in the form H = In <E> H± + H 2 ® In, 
with N x N matrices H\ and Hi- 

Matrices satisfy an infinite system of equations equivalent to the Eq. (1). These 
equations does not depend on the parameter A and can be consequently solved using the 
Buchberger algoritm ||. 

In [|j] this algoritm was succesfully applied to the solution of the constant Yang-Baxter 
equation with the .R-matrix does not depending on A. (Of course the condition (2) was not 
necessary in [|J). 

Our method of solving the Yang-Baxter equation is based on the well known fact that 
for all i?-matrices which were found up to now the quotents of their matrix elements may be 
rationally expressed on A, e A , or Jacoby elliptic functions on A. So studying the first terms of 
Taylor expansions for these quotents we may try to find them whole. 

The plan of the paper is the following. In the Sect. 2 we study the expansion (8) up to the 
order n — 6. In the Sect. 3 using the SL(2) symmetry (4) we separate all 4 x 4 Hamiltonian 
densities into standard classes. In the Sect. 4 using the equations presented in the Sect. 2 we 
find the complete list of integrable 4x4 Hamiltonian densities. In the Sect. 5 we obtain the 
corresponding i?-matrices. 

2 The Taylor expansion for R- matrix 

As it was mentoined in the previous section the first two terms of the Taylor expansion 
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are: 

R (o) = j R (i) = H (g) 

There are two ways to obtain the higher terms of the Taylor expansion of i?-matrix. First 
we may do it directly by substitution of partial sums into the Eq. (1). Another way is to 
differentiate the Eq. (1) by A and fi and then to take the limit A — * 0, /i — > 0. Technically 
the second method is more compact. However the first method gives a more detailed picture 
especially for n > 6 when in each order of A there are more than one equation. 

In the second order we obtain the following relation: 

itf - R$ = H\ 3 - E\ 2 . (10) 

Its solution will be = H 2 + 'jI. where 7 is an arbitrary numerical parameter. The term 7 J 
may be eliminated by the symmetry (3) if the function /(A) satisfy the following conditions: 
/'(0) = and /"(0) = —7. So we may write the following expression for R^: 

R^=H\ (11) 

In the order A 3 we obtain: 

R { S - R?2 = H 3 23 - H 3 12 + [(H 12 + H 23 ) , J] , (12) 
where J = [H 12 , H23]. The solvability of the Eq. (12) leads to the following condition: 

[(H 12 + H 23 ),J] = K 23 -K 12 , (13) 

where K is a A^ 2 x A^ 2 matrix. 

The equation (13) has the standard form: 

B = F 23 - F 12 , (14) 

with A^ 3 x A^ 3 matrix B and N 2 x A^ 2 matrix F The analogous condition also appears in the 
order A 5 . For studying the Eq. (14) let us consider F as an operator in C N (g> C N and B as 
an operator in C <g> C <8> C. For i — 1, 2, 3 let tr; be traces in the corresponding subspaces. For 
example if M x , M 2 and M 3 are some NxN matrices then tri (M 1 <g> M 2 <g> M 3 ) = trM x -M 2 (g)M 3 
and so on. 

If now for a given matrix B the Eq. (14) has a solution then the matrix F may be obtained 
by the following formula: 

F=±-tr 1 B + -^-tr 1 tr 2 B®I N + 1 I, (15) 
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where 7 is an arbitrary parameter. This result may be easely proved if we represent the matrix 
B in the general form: 

B = Xi <g> Yi + I N Zi + Ti <8> I N + aI N 2, (16) 

i 

where trAQ = trl^ = trZj = trTj = 0. 

Combining Eqs. (14) and (15) we obtain a system of linear equations on the elements of 
the matrix B. For the Eq. (13) this approach leads to the system of cubic equations on the 
elements of H . We shall solve it using the Buchberger algoritm. 

In the order A 3 we have met the first obstacle for integrability of the Hamiltonian density. 
If the matrix H depends on some parameters then solving the Eq. (13) we shall find the first 
system of restrictions on them. 

From (12) and (13) follows that: 

R (3) = h 3 + K. (17) 

The higher order expansion gives: 

= H A + 2(HK + KH), (18) 
R (5) = h 5 + L + 2(KH 2 + H 2 K) + QHKH, (19) 
R (6) = H 6 + KH 3 + H 3 K + 9(H 2 KH + HKH 2 ) + 

+ 10K 2 + 3(HI + LH). (20) 

The matrix L is defined by the formula: 

L 23 -L 12 = [H 3 2 + H 3 3 + 3(K 12 + K 23 ),J} + 3(H 12 [J,H 12 ]H 12 + H 23 [J,H 23 }H 23 ) + 
(H 12 H 23 + H 23 H 12 )(K 23 - K l2 ) + (K 23 - K 12 )(H 12 H 23 + H 23 H l2 ) - 
2(H 12 (K 23 - K 12 )H 23 + H 23 (K 23 - K 12 )H 12 ). (21) 

The formulas (19) and (21) follow from the relation: 

~ R-12 + 2(^12-^2? ~ ^12^23) + 3(^12-^2^ ~ H 23 ) = 

H 23 (R 12 R 23 )W - (R 12 R 23 )^H 12 , (22) 

which may be obtained from the Eq. (1) by differentiating it one time on fi and four times on 
A and taking the limit: A — > 0, fi — > 0. Here in (22) we denoted by (R± 2 R 23 ) ^ the following 
binomial sum: 

n 

( J Ri 2J R23) (n) = E^ jR 24r ) , (23) 
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where as usual ( = 7 . 5—7. 

n k\(n — ky. 

Using the similar procedure in the order A 6 we have also obtained the following equation: 

H 23 (R 12 R 23 )^ - (Ri2R2 3 ) {5) . (24) 
For even n the expressions for may be easely obtained from the equation: 

Rn(-V)RM = RwMRwi-v), (25) 

which follows from the Eq. (1) when A = 0. 

Let us notice that the Eqs. (20) and (24) are independent. So in the order A 6 the substi- 
tution of the Taylor sum for i?-matrix into the Eq. (1) gives two different equations. 

3 Classification of 4 x 4 Hamiltonian densities 

In this section we shall analyse the simpliest case of 4 x 4 matrices H. In each equivalence 
class under the action (4) of SL(2) group we shall choose a unique element. 
First let us represent an arbitrary 4x4 matrix H as the following sum: 

3 

H=^2 {SijOi <g> Oj + EijkAiGj ®a k + lif (I 2 <g> Oi + a t ® I 2 ) + (I 2 <g> a* - a t <g> I 2 )). (26) 
Here 5 is a symmetrical 3x3 matrix: 

S l = 5, (27) 

and A, and are vectors in C 3 . 

Matrices crj for % = 1...3 are the standard Pauli matrices: 

*-(;;)• *-(!?). -(i-0- <*> 

Using the data (S,A,h^ s \h^) and the well known fact of local isomorphism between 
5L(2, C) and 50(3, C) we may represent the SL(2) action (4) in the following equivalent 
form: 

5 = MSM~\ 
A = MA, 

U s > a) = Mh (s ' a) . (29) 
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Here M. is the corresponding to M 3 x 3 matrix satisfying the following equation: 

M l = M-\ 



(30) 



The Eqs. (27) and (30) may be represented in the following form: 



< Sx, y >=< x, Sy >, 



< Mx,My >=< x,y >, 



(31) 



where < ■, • > is a bilinear form on C 3 defined as: 
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< x,y >= ^xiyi. 



(32) 



i=i 



The form < •, • > is nondegenerative so that if for some ifC 3 and every y E C 3 < x,y >= 
then x — 0. 

Our study of the SL(2) symmetry will be based first of all on its action on the matrix S, 
whose standard form will be defined according to the structure of its eigenspace. Let k be 
the dimension of the space spanned on the eigenvectors of S. We shall separately consider 
the following three possibilities: k = 3, k = 2, and k = 1 and in each case define Sk as the 
standard form of the matrix S. 



In this case we have also to distinguish seven different subcases: 

1. All the eigenvalues are different and nonsero - the XYZ case. Then we shall use the 
following notations: V\ — v X i v 2 = v y , v 3 — v z . 

2. All the eigenvalues are different and one of them (v z ) is sero - the XY case. 

3. The two eigenvalues coincide [y x and v y ) and non of the eigenvalues is zero - the XXZ 



4. The two eigenvalues coincide and nonzero, however the third is zero [y x — v y — u, 
v z — 0) - the XX case. 

5. The two eigenvalues are opposite to each other and the third is zero [y x = v = —v y , 
u z = 0) - the X(-X) case. 

6. The two eigenvalues are zero and the third is nonzero [y x = u y = 0, v z = v ^ 0) - the Z 
case. 



3.1 k=3 



case. 



6 



7. All the eigenvalues coincide and nonzero - the XXX case. In this case v x — v y — v z — 
v± 0. 

8. All the eigenvalues are zero - the case. 

In all of this subcases we may find in the space C 3 the eigenbasis of the matrix S which 
will be orthogonal with the correspondence to the form < •, • >. When all eigenvalues of S are 
different this fact easely follows from the first Eq. (31) and nondegenerativity of < •, • >. When 
some of the eigenvalues are equal to each other we may use the ortogonalisation procedure. 
This also is possible due to the nondegenerativity of < •, • >. 

In the XYZ and XY cases the orthogonal eigenbasis of S is fixed. In this basis the matrix 
S has the following standard form: 



The Hamiltonian density depends on the nine parameters: z^, Ai, nf and h\ a . 

In XX Z, XX and Z cases there exist an additional symmetry group preserving the matrix 
S. It is generated by the following blockdiagonal matrices: 




v x 

Vy 

v z 




(33) 




(34) 



where the 2x2 matrix Mo satisfies the relation (30) and preserves the form: 



2 




(35) 



on C 2 . This means that for arbitrary x, y G C 2 . 



< M x,M y >o=< x,y > . 



(36) 



Using this symmetry we may reduce A acting on the 2- vector A = (A 1: A 2 ): 



A 







M A . 



(37) 



There are two different possibilities: 



1. < Aq, A >ot^ 0. Then A may be reduced to (Ai, 0, As). 



2. < Aq, A >= 0. Then A may be reduced to (A 1 ,iA i , A 3 ) 
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In the XXX case the symmetry group is the whole SL{2) and we also have to consider 
the two possibilities: 

1. <A,A >^ 0. In the appropriate basis A = (0, 0, A 3 ). 

2. < A, A >= 0. In the appropriate basis A = (A x , iA x , 0). 

3.2 k=2 

Let v\ and v 2 be the eigenvalues of S corresponding to eigenvectors e x and e 2 . We suppose 
that v\ is degenerative. This means that for some vector e$: 

Se = uie + e x . (38) 

We may find the triple (e , e\, e 2 ) satisfying the following system of relations: 

(ei,e fc ) = 0, k = 1,2, (39) 
(e ,e 2 ) = 0. (40) 

When vi ^ v 2 the Eqs. (39) and (40) are satisfied automatically. The condition (e 1 ,e 2 ) = 
follows from (31). Also 

(ei, ei) = ((S - v x )e Q , e x ) = (e , (S - u^d) = 0, (41) 

and 

(ei, e 2 ) = ((5 - i/i)e , e 2 ) = (f 2 - fi)(e , e 2 ) = 0. (42) 

For V\ = v 2 the condition (40) is more refined. From the Eq. (39) and nondegenarativity 
of < •, ■ > follows that (e , e x ) ^ 0. So if the vector e 2 does not satisfy the condition (40) we 
may define a new one 

(e ,e 2 ) 

e 2 = e 2 -- r-ei. (43) 

(e ,ei) 

Now all the relations (39), (40) are satisfied for the triple (e ,ei,e 2 ). 

The 3x3 symmetrical matrix 5 whose eigenvectors satisfy the conditions (39) and (40) 
has the following canonical representation: 



/ fx + f -5- U \ 

-5- v x - u 
V ' v z / 



(44) 



where z/ x = 1/1, z/, = z/ 2 and Q = ; 6 °' ^ . 

* (eo,eo) 

As in the previous case we need to consider now five different subcases: 
I. The XXZ case - eigenvalues are nonsero and v x ^ v z . 
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2. The XX case - u x = v ^ 0, v z = 0. 

3. The Z case - v x — 0, v z — v ^ 0. 

4. The XXX case - i/ x = i/ z = v. 



5. The case - v x — v z — 0. 



3.3 k=l 

First let us consider the following matrix 

S — S — ul^, 



(45) 



where v is the eiqenvalue of S. 

The matrix S is double degenerated. This means that there exist a triple (e_i,e ,ei) with 
the following properties: 



Se\ = 0, S^o = ei, .Se-i = eo- 



(46) 



The matrix S 2 has two different eigenvectors e\ and eo with zero eigenvalues. According to 
(44) it may be reduced to the following form: 



1 % 
S 2 = f3 2 \ i -1 


where (3 is a new parameter. The vectors e±, eo and e_i may be represented as follows: 

/ 1 \ 


5 ' 

v J 

where p and 6 are new parameters. 

The system (27), (46) defines the explicit form of S: 

5 i5 (3 
S = I i5 -5 i(3 

P iP 



(47) 













ei = /3 2 








e-i 




vv 









(48) 



(49) 



as well as the condition p — p. Now the matrix S may be represented to the following canonical 
form: 

v + 5 i5 p 
S s = I i5 is -6 ip ] . (50) 

/3 ip v 
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As in the previous cases we shall consider the two different subcases: 
1. The XXX case v ^ 0. 



2. The case v = 0. 

4 Classification of integrable Hamiltonian densities 

In this section we present the complete list of 18 classes of integrable 4x4 Hamiltonian 
densities. It is a remarkable fact that for obtaining all of them it was necessary to use only 
the Eq. (13). The Eqs. (21) and (24) have satisfied automatically. We may contend that 
these solutions really correspond to integrable matrices H becouse in the next section the 
corresponding /^-matrices will be presented. 

The result was obtained by computer calculation using the Buchberger algoritm. We have 
used for this purpose the computer algebra system MAPLE V, Release 3 ||,||. The result 
will be represented as a list of the following sets: G = {gi, g n }. Each of them is equivalent 
the solution: gi — ... — g n — 0. In our calculations we have used the command gsolve() of 
the Grobner package. In order to present solutions in more compact form we have introduced 
the following variables: h±' a ^ = h^'^ ± ih^ , A± — A\ ± iA 2 . Sometimes we have used the 
following notations: h^i'^ = /ij and A z = A3. Also we have used the standard shortening 
changing for example A\\ A 2 ; A 3 by ^1,2,3 • 

As in the previous Section we shall separately consider the three possibilities: k = 3, 2, 1. 



4.1 



k=3 



4.1.1 



The XYZ case 



£3,1 - {^1,2,3; ^2,3; ^1?2,3}- 



(51) 



4.1.2 



The XY case 



(^3,2 — {A, 2,3; h^; h± 2 3 }. 



(52) 



4.1.3 



The XXZ case 



£3,3 ~~ {Ai.,2; ^1*2,3! ^1,2}- 



(53) 



4.1.4 



The XX case 



^3,4 _ {^1,2; ^1*2) ^1?2,3}- 



(54) 
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4.1.5 The X(-X) case 

4.1.6 The Z case 

Exept the degenerative case of 6*3,3 there are two new solutions: 

G*3,6 — {^-1,2,3; h^; ^^2,3} i 
G 3 j = {^1,2,3; h^ 2 ; h^}. 

4.1.7 The XXX case 

In this case all the solutions are the degeneratine cases of 6*3,3 an d G 2jG . 

4.1.8 The case 

Exept the degenerative case of G 2 j there is one solution: 

4.2 k=2 

4.2.1 The XXZ case 

. /» ■ /» 



G 2 ,l - {^1,2,35 ^1*2,3! ^1?2,3) 



4.2.2 The XX case 



G 2 , 2± = {A 1 , 2 ;h[%h^v±hi s) -2A 3 ±ih<? ) }, 
^2,3 = {^1,2,3; M^! ^1^2,3}- 



4.2.3 The Z case 



G 2A = {A + , 3 ; h%; h{% av + (fc«) 2 ; ah? - iA.hP) 
G 2 ,5 = {^1,2,3; ^+,3; ^"2,3}- 



4.2.4 The XXX case 



G 2 ,6 — {-4+,3; ^1,2,3! ^+,3}- 
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4.2.5 The case 



G 



2,8 — 



{A + ;h[% 3 ;h^;h^A 3 -h^A_}. 



(65) 
(66) 



4.3 k=l 



4.3.1 The XXX case 



- {^1,2,3! ft%2,3; ^2,3} • 



(67) 



4.3.2 The case 



Gi, 2 = {A +i3 ; h%; h^;8h 3 a) + (3h w - iA_h { l ] }. 



(68) 



5 The corresponding i?-matrices 

In this sections we present i?-matrices corresponding to all the solutions presented in the 
previous sections. For some of the i?-matrices the parameter A is as vector: A = (Ai, A2, A n ). 
This happens when the dependence of elements of the i?-matrix from A has the following form: 



If all the parameters Oj are independent we may put A; = a^A and consider A as a vector 
parameter. 

According to their dependence on A all the presented i?-matrices may be separated in 
four classes: 1. elliptic, 2. trigonometric, 3. rational-trigonometric, and 4. rational. When 
to different values of parameters of Hamiltonian density correspond i?-matrices of different 
classes we have used the additional labeling. The trigonometric /^-matrices were labeled by 

the rational-trigonometric by Hr\ and the rational by V. 



Ri,j(X) = RijififaX), / 2 (a 2 A), f n (a>nX))- 



(69) 



5.1 k=3 



5.1.1 The XYZ case 



The i?-matrix corresponding to G31 is the Baxter elliptic solution of the Eq. (1) 
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5.1.2 The XY case 



The i?-matrix corresponding to G3 2 is the Felderhof elliptic solution of the Eq. (1) M. 
special case it degenerates into the following trigonometric f?-matrix: 



R 



3,2t 



( shr] chA + shA shr? shA \ 

shr/chA ch?7shA 

chr/shA shr^chA 

V shr/ shA sh?7 chA — shA J 



(/»W) a 



v x v v . 



5.1.3 The XXZ case 



R. 



3,3i 



R 



3,3tr 



( sh(Ai + rj) \ 

e A2 shr/ 6»shA! 

y shAi e" A2 shr] 

sh(Ai + 7?) / 








V 

/ Ai + 77 

rje X2 0Ai 

^ r/e" A2 

y 



\ 






5.1.4 The XX case 



R. 



3,4t 



3,4tr 



/ sh(Ai + 77) 

e A2 shr/ 6>shAi 

^shAi e A2 slrr/ 

V 

/ A1 + 77 






V 







r/e A2 6>Ai 
^ ^e A2 









-A1 + 77 / 



-sh(Ai - 77) / 



5.1.5 The X(-X) case 



R 



3.5 



/ chA x sinA 2 \ 

C0SA2 shAi 

— shAi COSA2 

y sinA 2 chAi J 
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5.1.6 The Z case 

For htphP ^ the solution G 3j6 can be reduced by the SX(2)-symmetry to the Felderhof 



elliptic case. To the case h+ h_ = correspond matrices i?3,6± (h±' = 0). 



(s) 



R 



3,6+ 



#3,7 



/ 1 9(1 - 


e A ) 


9(1 - 


e A ) 2#(chA - 1) \ 


e A 









9(1 - e A ) 







e A 




9(1 - e A ) 


\ 









1 / 


/ e A i 








\ 




e A2 















e A 3 









\ 





e A 4 


/ 





-^3,6+' 



(76) 



(77) 



5.1.7 The case 



R. 



3,8t± 



R 



■3,8tr± 





sh(Ai + 77) 















sh?7e A2 


=FzshAi 












±ishAi 


sh?7e~ A2 







V 













sh(Ai 


/ 


Ai + V 

















r]e X2 =rAi 














±iAi ?7e~ A2 







3 


\ 








-Xi + 


V / 





W\2 
3 ) 



(7f 



(79) 



5.2 k=2 
5.2.1 The XXZ case 

/ sh(A + 77) sh(A + 7/) shA \ 

shr? shA 

shA sh?7 

V sh(A + r/) y 

In the special case: 9 = — 4shr7 this i?-matrix was obtained in [|J. 



R2A 



(80) 



5.2.2 The XX case 



R 



2.2t± 



( sh(A + 77) 

sh?7e~ A ±e _r? shA 

ie^ shA shr] e A 

\ 



6>shA \ 



-sh(A — 77) J 



14 



R 



2,2r± 



R 



2.3 





V 

±A 




/ A + r/ 




/ sh(A + r/) 



V o 



o xe \ 
±x o 

7] o 

-A + 7/ / 

#sh2A 

shi] shA 

shA shr] 

-sh(A - 



A* = h 



\ 



V) J 



(a) 



5.3 The Z case 



ft 



2,4t 



2Ar 



R 



2,5 



/ e Al #(e Al - e" A2 



A 2 



Ai 

e" A2 



\ 

/ 1 ex -9X X(6 2 X + 1) \ 

10 9X 

001 -ex 

\ 1 

/ „A a „r\ r ,u\ a „w r ,u \ e 2 



6»(e Al -e A2 ) 2e 2 (chX 2 




0(e Al - e" A2 ; 
#(e Al - e Aa ) 



/ 



,(«) 



0. 



e A fle^shA tfe^shA ^- (ch2r/ e 3A + sh2r/ e" A 






V 



5.3.1 The XXX case 












0e"shA 
0e"shA 



R 



2,6 



X + T] 




-6 A 


#A + ^+^A 2 







A 


6A 





A 


'/ 


-e 2 A 











A + r/ 



/ 

In the special case: £i = £2, # = -£162 this .R-matrix was presented in 



5.3.2 The case 



/ 1 


(^ + 0(e A 


-1) (6 




-e" A ) AshA + 2(fl 2 + £r/)(ch\ 





e A 









(e + V )(e x -l) 










e" A 




(e-7])(l -e" A ) 


\o 












1 


/ 1 


(0 + OA 


[e-ox 


AA + ( 


9 2 + ir ] )X 2 ) \ 







1 





(6 


+ T))X 


, /4 a) = 0, 








1 


(6 


-T))X 












1 ) 
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/ chAi 6>(e A2 -e- Al ) ^(e" A2 -e Al ) r^shA 2 + 2# 2 (chAi - chA 2 ) \ 



R 



2.8 






V o 



e A 2 

shAi 




-shAi 

e -A 2 





0( e Ai _ e A 2 ^ 

0(e 



-Ai _ e -A 2 -j 

chAi 



5.4 k=l 
5.4.1 The XXX case 

( X + rj 9X(X + rj) 9\(X + ri) -^A(A + V )(X 2 + A?? + £) ^ 

77 

A 

\ A + 77 / 



#1, 



A 

'/ 




-9X{X + V ) 
-9X(X + V ) 
A + 77 



In the special case 9 = — 1, i = rf this _R-matrix was obtained in ||. 
5.4.2 The case 

In this case we have separately considered the two subcases: 8h^ 7^ and 8h% 
first subcase: 



Ri, 



( 1 ^i(l -e" A ) 


^ 2 (e A - 1) 


2# 2 (A(chA - 1) - 


shA) \ 


e" A 





^3(1 -e" A ) 







e A 


^ 4 (e A - 1) 




\ 
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/ 



where 



In the second subcase: 
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e Vl eh?? + e' 
e^chr/ + e _</32 
e^dxrj - e~^ 2 
'sht] - e _</31 
sh((^i + c^ 2 )sh2^. 



pV2 



! sh?7, 
! sh77, 
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e" A 
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|(e" A -l) 



5 = 0, 



R\,2r — 



( 1 0f A 9r]X -9 2 X 2 + 7 A \ 
1 








1 





-|a 

1 



hf ] = 0. 



(95) 



/ 



Strictly speaking the list of presented i?-matrices is not full. For example making the 



contraction: 9 — > oo, rj 



7, we may obtain from R^^t the following i?-matrix: 



R 



(96) 



/ e Al \ 

e A2 7e Al 
3,3* - o e~ A2 
\ 

This case corresponds to u x = %A Z and it also may be easely be obtained from the following 



e Al / 



expressions for 9 and r\: 9 



6 Conclusion 



Z/7. ~\~ 1 

- x — ■ — ^-r 2 -, shr] 
v x — IA Z 



In this paper we have presented a " direct" computer algebra algoritm for solving the quantum 
Yang-Baxter equation using the Taylor expansion of i?-matrix. Applying this algoritm in the 
4x4 case we have obtained a complete list of solutions. 

We beleive that our method will be useful for studying of concrete Hamiltonian densities 
appearing in various branches of theoretical physics. 



References 

[1] P. P. Kulish, E. K. Sklyanin 1980 Solutions of the Yang-Baxter Equation Zapiski Nauch. 
Sem. LOMI95 129 [1982 J. Sov. Math. 19 1596] 

[2] V. E. Korepin, A. G. Izergin, N. M. Bogoliubov 1993 Quantum inverse scattering method 
and correlation functions (Cambridge Univ. Press) 

[3] B. Buchberger 1982, 1983 Grobner Bases: An Algoritmic Method in Polinomial Ideal 
Theory in Computer Algebra Symbolic and Algebraic Computations (Springer) 

[4] J. Hietarinta 1993 Solving the two-dimensional constant quantum Yang-Baxter equation 
Journ. Math. Phys. 34, 1725 

[5] K. M. Heal, M. L. Hansen, K. M. Rickard 1996 Maple V Learning Guide (Springer) 



17 



[6] M. B. Monagan, K. O. Geddes, G. Labahn, S. Vorkoetter 1996 Maple V Programming 
Guide (Springer) 

[7] B. U. Felderhof 1973 Diagonalization of the transfer matrix of the free-fermion model II 
Physica 66 279 

[8] I. V. Cherednik 1980 On a method of constructing factorized S-matrices in elementary 
functions Teoret. Mat. Fiz. 43 li7[1980 Theoret. Math. Phys. 43 356] 



[9] P. P. Kulish and A. A. Stolin 1997 Deformed Yangians and integrable models q-alg/9708024 



18 



